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LONG WRITE-UP 



1. Introduction 

The success of current and future collider experiments depends on a pre- 
cise prediction of the expected cross-sections for both signal processes and 
their corresponding backgrounds. The experimental precision of these col- 
liders can only be matched by theoretical calculations beyond leading order 
in perturbation theory. Amongst the most pressing problems are the vir- 
tual (one-loop) corrections to processes with up to four particles in the final 
state. A disection of these one-loop amplitudes into gauge invariant sub- 
amplitudes and the use of a compact notation allows for the generation of 
numerically stable, fast computer programs for their evaluation. The authors 
of [H introduced the spinor helicity formalism for massless particles in their 
cross-section calculations. Together with its extensions for massive particles, 
this approach leads to a compact representation of helicity amplitudes. 

Although recently a number of promising purely numerical methods for 
the calculation of one-loop amplitudes have been presented 



semi-nu- 



merical and algebraic methods are still an important tool for matrix element 



calculations |17H37I|. The interested reader finds a more complete list of 



methods for one-loop ampitudes and their applications in [38|, |39 . 

The algebraic manipulations required in matrix element calculations are 
very often simple, local operations that have to be applied to millions of 
terms, not requiring the knowledge of the whole expressions at any given 
moment^. General purpose computer algebra programs easily hit the memory 
limits when dealing with expressions of this size as they hold the whole 
expression in memory at every point in the program. A different approach 
is persued by the symbolic manipulation program Form [40]. The program 
provides only local manipulations of expressions, a fact which is reflected 
by the memory model, viewing an expression merely as a stream of terms; 
only a single term needs to reside in memory at any moment while the rest 
of the expression is stored on a storage device. In practise, a system of 
buffers is used to reduce the number of disk operations and multiple terms 
are processed at once on multi-processor systems. 



*An example for a non-local algorithm is the factorization of a polynomial, which needs 
to know all its terms at once. 
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The spinor helicity formalism has been implemented in the computer 



algebra system Mathematica through the package S@M |41[. It provides 
routines for the algebraic manipulation and numerical evaluation of spinor 
products and Dirac matrices. Mathematica's numerical and algebraic capa- 
bilities paired with a user-friendly interface are a clear advantage for many 
calculations. However, the above mentioned memory limitations restrict its 
usage to smaller problems. More involved cross-section calculations require 
an implementation of spinors in Form, combining the elegance of the spinor 
helicity formalism with Form's ability to process huge expressions. In this ar- 
ticle we demonstrate that it is very easy to exploit Form's existing constructs 
to implement helicity spinors. Therefore the implementation is in the guise of 
a Form library, spinney, rather than an extension of the core language. The 
naming of many of the functions and procedures makes reference to S@M 
allowing the user an easier migration between the two libraries. 

This article is structured as follows: Section [2] gives an overview over the 
underlying theory and establishes the notation used in the rest of this work. 
Section [3] provides an interface documentation of the provided procedures 
and the defined symbols and functions. A couple of examples are discussed 
in Section HI 



2. Theoretical Background 

2.1. Conventions 

We consider spinors u(p) and v(p) which are solutions of the Dirac equa- 
tions 

(f — m ■ T)u(p) = and {f + m ■ I)v (p) = 0. (1) 

Here, we denote I the identity operator in spinor space and p 1 = g^^p" ■ 
Dimension splitting is understood in the way that for n e C 

ST = r + ST (2a) 



such that 



n — 4 and 



r P ~9 P u = 0. 



(2b) 
(2c) 
(2d) 
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The projector into the four- dimensional, physical Minkowski space is of the 
form 

$ = diag(l,-l,-l,-l)©0< n - 4 \ (3) 

i.e. it leads to the on-shell condition p 2 — m 2 = 0. 

As a short-hand notation we add twiddles and hats to all vector-like ob- 
jects to indicate projections into the 4 and (n — 4) dimensional subspaces 
(7 M = g^ u , = g^k u , etc.). We work with anti-commuting 75 in 4 dimen- 
sions, which means 

{7", 7"} = 20"" and (4a) 

{75, r} = [75,71 = 0. (4b) 

The characteristic equation (7s) 2 — 1 = allows to introduce the usual pro- 
jectors and Ur into the left-handed and right-handed subspace, 

n L = n_ = i(I- 75 ) and U R = U + = ^(I + l5 ). (5) 

The helicity eigenstates of massless spinors are denoted by a commonly used 
bracket notation^ 

U + u(ki) = H + v(ki) = \i), U^u(ki) = n„t>(A;j) = (6a) 

u{ki)U_ = v(ki)IL- = [i\, u{ki)U + = v(ki)U + = (6b) 

Also in the massive case we use a similar bracket notation to distinguish 
the different solutions of the Dirac equation. We notice that using a given 
lightlike vector q every massive vector pj can be decomposed into a sum of 
two lightlike vectors as in the following equation 

which defines the lightlike vector k{. The solutions of the Dirac equations 
{jj>i ± m/)!/ 1 * 1 ) = and {j>i ± m/)]/*] = can be expressed in terms of the 



2 We identify the vectors with their labels where this does not lead to ambiguities, e.g. 
pi = I or kj = j. Throughout this paper lower case Latin labels are used for light-like 
vectors and upper case Latin labels for massive vectors. 
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massless spinors of q and fcj, 



|/±> = |i>± M 1 " 1 ' |/±] = {i] ± ft 1 *' (8a) 

<i±l = «l±^bl, M = W±§<«I- < 8b > 

It should be noted that in spite of the similarity in notation the massive 
spinors are not constructed as eigenstates of the helicity projectors U±. 

2.2. Dimension Splitting 

In this section we want to derive the formula? required for separating the 
dependence on g from the spinor chains and spinor traces. The starting point 
for such a separation is the following equation, a proof of which can be found 
in 
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tr{I} tr{7 w • • • 7^7^ • • • 7^} = tr{ 7 w • • • 7^} ti{^ ■ ■ ■ 7^} (9) 

A chain of Dirac matrices can always be sorted by splitting 7^ = 7^ 4- 7^ 
and shuffling 7 and 7 to opposite ends of the chain using Equation (JH). The 
presence of 75 does not change the above results as in our scheme it can be 
written as 75 = ie PLVpa Yl u l p l a '/4!. 

Equation (Q extends to the case of spinor chains by the observation that 
for any pair of massless spinors delimiting a product of Dirac matrices one 
can introduce light-like auxiliary vectors p and q in order to turn the spinor 
chain into a trace. 



_ tr{n_fc---^} ... , . tr{II_fc---fcM 

\J\ — lAU\i\ ' \ l \"'\J) — 



UWM 

tr{n + fc---^} r| , , tr{n + ^---^M 



(10) 



In all four cases, the insertion can be undone after Equation (Q has been 
applied and hence no spurious denominators need to be inserted in an actual 
calculation. 

The trace of (n — 4) dimensional Dirac matrices 7 can be evaluated ac- 
cording to the recursion relation 

tr {f 1 . ..y*} = J2 = 2 q (-l)Y lVi tr{f 2 • • .7^-17^+1 ■ ■ .7"*} . (11) 
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Since only even numbers of 7 lead to non-vanishing contributions, a chain 
with two Dirac matrices can be expanded according to the above rules into 

= (i\rr\j) + wrrij} = + r v m (12) 

As a more complicated example, we also give an explicit expression for a 
spinor chain with four matrices: 

+ r p [Arr\j} + «r wrr bl - sr wrr bl + sr wrr bl 

+ [u] (fiTsT - S^fiT + sTsT') • (13) 

,2.5. Chisholm Identities 

With the same argument as in the previous section, we can extend the 
validity of the Chisholm identity 

tT{Y L ■ ■■T lm ~ x l v }% = 2 (7 W • ■■Y 2m ~ 1 +Y 2m ~ 1 ■ ■ - 7 m ) (14) 

to the case of spinor chains by substituting Equation (JT0|) . If we denote 
products of Dirac matrices as T = 7 m • • • 7^ and T' = 7" 1 • • • j"* and their 
reversed strings by r 7 = 7^ p • • • 7 /il , the Chisholm identities for spinor chains 
read 



(i\Vrr'\j) ■ % = 2V\j)(i\T - 2f\i)(j\V (15a) 

[z|r 7 M r'|i] • % = 2r'\mr - 2r T K][j|r 7 (isb) 

(i\Vrr'\j} ■ % = 2r'\j]{i\T + 2 f\i)[j\T' (15c) 

[i|rfT / |j> • % = 2r'| J )[z|r + 2i T K](j|r 7 (i5d) 



The above identities are valid only if the number of Dirac matrices in T and 
T' matches the helicities of the spinors % and j, i.e. in equations fll5ap and 
( 115bj) the length of r^T' must be everJl and in (I15cp and ( 115dj) it must be 



odd. 

The repeated application of Equation ffl5|) together with the identities 



j»T% = -2f and f t Tr% = 2(YT + fr) (16) 



75 counts as an even number of matrices 
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for an odd number of matrices in T, ensures that, starting from an expression 
with no uncontracted Lorentz indices, all spinor chains in four dimensions 
and traces can be expressed in terms of spinor products of the form (ij) 
and [ij]. 

2.4- Majorana Spinor s 

In supersymmetric extensions of the standard model we need to deal 
with interactions of Majorana fermions. A Majorana fermion is its own anti- 
particle, i.e. it is invariant under charge conjugation: 

■0Af = C^Pm = (17) 

The charge-conjugation matrix has the properties: 

Ct = C'\ C T = -C, CTfC" 1 = rjiTi (18) 

with 

We aim to write down a consistent set of Feynman rules to deal with Ma- 
jorana fermions. The problem is that vertices involving Majorana fermions 
violate fermion number flow; the fermion flow in the Feynman diagram is ill- 
defined. A consistent way of dealing with this was proposed in j43[ . There are 
two drawbacks with this approach; charge conjugation matrices are explic- 
itly introduced into the Feynman rules and the relative sign of the Feynman 
graphs needs to be determined from the original Wick contractions. 



We follow the approach in |44j which has been implemented in (45 
Each vertex containing Dirac fermions has two expressions, one in which the 
fermion flow follows the fermion number flow, and the other "flipped" ver- 
tex, where the fermion flows in the opposite direction to the fermion number 
flow. In our implementation we impose a fermion flow on the vertex through 
the procedure RemoveNCContainer and "flip" the vertices and spinors where 
neccesary. 

We have implemented the following flipping rules, using f fT9"j) : 

W ' (20) 

(n ± )' = n ± (20) 
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and for the spinors 

\P)' = (P\ (<y,\ 
\p)' =[P\. [ ] 

This leads to the following rule for the fermionic propagator: 

S'(p) = — L_ = S(-p). (22) 
—p — m 

Another appealing feature of this method is that the relative sign be- 
tween graphs can be determined directly from the expressions as opposed to 
reverting back to the original Wick contractions. This brings it in line with 
the usual formulation for Dirac fermions. To compute a consistent relative 
sign of Feynman graphs each must be multiplied by (— 1) P+L , with, 

• P: the parity of the permutation of external spinors with respect to 
some reference order, 

• L: the number of closed fermion loops. 

This sign is determined after the flipping rules are applied. 

3. Program Description 

3.1. Installation 



The library has been written using the literate programming tool nuweb 46 
It can be obtained by downloading the file spinney. tgz from the URL 
http : / /www . nikhef . nl/~thomasr/ f iletransf er . php=spinney . tgz . The 
tarball contains the following files: 

spinney. hh the Form file 
spinney.pdf annotated source code 
spinney. nw the nuweb sources 
spinney_test . f rm unit tests 

After unpacking the library the should run the test program spinney_test . f rm 
in order to ensure that the installed Form version is recent enough. If one or 
more tests in the program fail the most likely reason is a deprecated version 
of Form. 

The file spinney, nw is only needed if the user wants to rebuild any of the 
other files from scratch. 
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3.2. Representation in Form 

The spinney library uses four different representations for spinorial ob- 
jects, each of which is considered the most convenient form at a given point 
in the workflow. 

For the description of the definitions and declarations we use the Form 
keywords (CFunction, NFunction, CTensor, NTensor, Vector, Sym- 
bol, ^Procedure) in the text to indicate its type. Where more than one 
definition is combined in one box the arguments apply to all of the definitions. 

3.2.1. N on- Commuting Objects 

The 't Hooft algebra and the dimension splitting are applied to Form's 
non-commuting object. In order to keep the information about the dimension 
in the function names rather than the indices we do not use Form's Dirac 
matrices here (g_) but introduce our own functions. 

NTensor Sm(/i) 

An n-dimensional Dirac matrix 7 M . 
Parameters 

/i : Lorentz index or vector 

NTensor Sm4, SmEps(/i) 

The 4-dimensional and (n — 4)-dimensional (resp.) projections 
7 M and 7^ of the Dirac matrix 7^. 

The non-commuting objects corresponding to the four-dimensional ma- 
trix 75 and the derived projectors TT_|_ = (I±7 5 )/2 are described below; they 
have been implemented as non-commuting functions without parameters. 

NFunction Gamma5() 

The Dirac matrix 75. 

NFunction ProjPlus, ProjMinusQ 
The projectors Ur = n + = (I + 75) /2 and 

II L = II_ = (I - 75 )/2 (resp.) 

In the case of the spinors we use the same function names for both mas- 
sive and massless spinors. Massive spinors are distinguished by their second 
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argument p, which is ±1, indicating that the spinor lies in the kernel of 
the operator {j> — pm). For massless spinors the second argument must be 
omitted. 



NFunction USpa, UbarSpb(p [, p]) 

The massless spinor \p) and its conjugated [p\, if p is omitted; 
the massive spinor \p p ) and its conjugated [p p \ otherwise. 

Parameters 

p : momentum of the spinor 

p: (optional), a value ±1 as described above 

NFunction USpb, UbarSpa(p [, p\) 

The massless spinor \p] and its conjugated (p\, if p is omitted; 
the massive spinor \p p ] and its conjugated (p p \ otherwise. 

Parameters 

p : momentum of the spinor 

p: (optional), a value ±1 as described above 

The following two functions mark the beginning and the end of a Dirac 
trace. 

NFunction trL, trR() 

Indicates the begin (resp. end) of a Dirac trace. 

For example, the expression tr{II + 7 M ^7^} would correspond to the fol- 
lowing product in Form: 

trL * ProjPlus * Sm(mu) * Sm(p) * Sm(nu) * Sm(q) * trR 

3.2.2. Indexed Notation 

In this form the spinor indices of Dirac matrices and spinors are kept. 
For example, the term (/cilp^A^] would be expressed as 

a 2 

and in a Form program as 
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• \ 

NCContainer (UbarSpa (kl ) , alphal) * 

NCContainer (USpb (k2 ) , alpha2) * 

NCContainer ( Sm (p ) , alphal, alpha2 ) 

V J 

This notation is particularly useful when importing expressions from a 
diagram generator, which does not necessarily put the factors in the cor- 
rect order. All elements of the spinor line are wrapped inside the function 
NCContainer which is defined as follows. 

CFunction NCContainer (o, %\ [, i 2 \) 

Representation of a spinor or a Dirac matrix with explicit 
spinor indices. 

Parameters 

o: a non-commuting object, such as Sm, Gamma5 or 
UbarSpa; a complete list is given in 
Section 13.2.11 

ii : for spinors: the only spinor index; for Dirac 
matrices: the first spinor index 

i 2 : for spinors: not present; for Dirac matrices: the 
second spinor index 

3.2.3. Collected Form 

The collected form is generated from non-commutative objects by the pro- 
cedure SpCollect or from the open form by the procedure SpClose. Spinor 
strings are represented by the following four functions 

CFunction Spaa, Spab, Spba, Spbb(/c 1? . . . , fc 2 ) 

A string delimited by two spinors. The four forms correspond 
to . . . \k 2 ), (ki\ . . . \k 2 ], [h\ . . . \k 2 ) and [k x \ . . . \k 2 ] 
respectively. 

Parameters 

ki : a light-like, four-dimensional vector. 

. . . : a list of four-dimensional indices or vectors. 

k 2 : a light-like, four-dimensional vector. 
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3.2.4- Open Form 

The open form can be generated from the collected form by calling the 
procedure SpOpen after all Lorentz indices have been contracted using SpContract. 
The reverse operation of SpOpen is SpClose, which transforms the open form 
back to collected form. 

CFunction Spa2, Spb2(&i, k 2 ) 

The spinor products (kik 2 ) and [kik 2 ] respectively. These 
functions are defined as anti-symmetric in their arguments. 

Parameters 

ki : a light-like, four-dimensional vector. 
k 2 : a light-like, four-dimensional vector. 

3.2.5. Metric Tensors 

The dimension splitting of g^ u into the sum g^ + as described in 
Equation fl2]) is reflected in the definition of the metric tensor and its projec- 
tions onto subspaces. As in the case of the Dirac matrices we find it more 
convenient to define new functions for the metric tensor rather than using 
the tensor d_, which is predefined in Form. 

CTensor d(/x, u) 

The n-dimensional metric tensor g^ v . This function is defined 
symmetric in its arguments. 

CTensor d4, dEps(/i, v) 

d4 is the 4-dimensional projection g^ v of the metric tensor; 
dEps is the orthogonal projection g^ u into the (n — 4) 
dimensional subspace. 

These function are defined symmetric in their arguments. 

3.2.6. Other Symbols and Objects 

The function SpDenominator has been introduced in order to represent 
the reciprocal value of its argument. Therefore, the replacement SpDenominator(x 
1/x is always safe. However, keeping certain denominators inside function 
arguments makes some substitutions easier. 
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CFunction SpDenominator(x) 

Represents the reciprocal value of the argument, i.e. 1/x. 

The function SpERRORTOKEN has been introduced to indicate inconsiten- 
cies detected by the procedure SpCheck which signals if the length of a spinor 
string trivially nullifies an expression. Although these null-expressions are 
not errors, in certain circumstances they can indicate errors in a program. 

CFunction SpERRORTOKEN() 

Indicates that the procedure SpCheck found an error in the 
term where SpERRORTOKEN appears. 

3.3. Reserved Symbols 

The library spinney defines more objects for internal use. All imple- 
mented algorithms require that none of these objects are present in any 
active expression at the invocation of the procedure and ensure that these 
objects are not present in any active expression after the procedure returns. 
These objects are declared exactly as written below. 



CFunctions fDUMMYl, 


. . . , fDUMMY4; 


Symbols sDUMMYl, . . 


, sDUMMY4; 


Indices iDUMMYl, . . 


, iDUMMY4; 


Vectors vDUMMYl, . . 


, vDUMMY4; 


NFunctions nDUMMYl, 


SpFlip ; 



We also define the three sets SpORIGSet, SpIMAGSet and SpObject which 
are used inside RemoveNCContainer. 



3.4- Implemented Algorithms 
3.4.I. Light- Cone Decomposition 

This algorithm uses Equations (JTj) and (IE]) in order to express massive 
vectors and spinors in terms of massless ones. 
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^Procedure LightConeDecomposition(jor, pi, q, mi) 

In every active expression all appearances of pi and all spinors 
|J ± ), | J ], (J ± | and [7 ± | are replaced by a pair of light-like 
vectors (resp. spinors). 

The parameters are defined according to Equation ([7]): 
Parameters 

pi: a massive vector that fulfills (pi) 2 = m 2 - 

Pi : a light-like vector 

q : a light-like vector (reference momentum) 

mi : the mass of pi 



This procedure works on non-commuting objects. If massive 
spinors are present this procedure must be called before 
SpCollect. 



This procedure needs to be called before SpCollect as in the closed form 
only massless spinors are allowed. 

As an example, we write down the numerator of the color stripped tree 
level diagram of uu — > ti in QCD for two different helicities. 



Local dl 
Local d2 



UbarSpa(kl 
UbarSpa (p3 
UbarSpajkl 
UbarSpb (p3 



* Sm4(mu) * USpb(k2) * 

+ 1) * Sm4(mu) * USpa(p4, 

* Sm4(mu) * USpb(k2) * 

+ 1) * Sm4(mu) * USpa(p4, 



■i); 



T^call LightConeDecomposition (p3 , k3 , kl , mT) 

T^call LightConeDecomposition (p4 , k4 , k2 , mT) 

#call SpCollect 

#call SpContract 



Print ; 
. end 

We have also added the command SpContract in this example as it simplifies 

the output considerably, which is given below. 

i >i 



15 



dl = 2*Spaa(kl ,k4)*Spbb(kl ,k2)* 

SpDenominator ( Spb2 (kl , k3))*mT; 
d2 = 2*Spaa(kl,k4)*Spbb(k3,k2); 
0.00 sec out of 0.03 sec 



The above calculation can therefore be summarized as: 



di = (£:i|7 M N(P3 \%M 
3.4-2. n- Dimensional 't Hooft Algebra 



(kik 4 )[kik 2 ] 
2 ■ r . . ; ■ rriT 



[kih] 

{W*\ki](pt\%\Pl] = 2 • {k x k A )[k z k 2 



^Procedure tHooftAlgebra() 



Carries out the (n — 4)-dimensional part of the algebra such 
that only four dimensional Dirac matrices 7 M are left; the 
dependence on (n — 4) is entirely expressed in terms of g^ u . 

This procedure works on non-commuting objects. If 
n-dimensional Dirac matrices, or 75 are present this 
procedure must be called before SpCollect. 



As an example, we give a short program that reproduces Equations ((T2 
and (H5D. 



Local examplel 

UbarSpa(kl) 
Local example2 = 

UbarSpa(kl) * Sm(mu) * 
Sm(rho) * Sm(sigma) 
#call tHoof tAlgebra () 
#call SpCollect 
Print +s ; 



Sm(mu) * Sm(nu) * USpa(k2 



Sm(nu) * 
* USpa(k2 



3-4-3. Change of Representation 

As described in Section I3.2[ different representations for the expressions 
are used at different points in a Form program using spinney. Here, we give 
an overview of the functions changing between the representations. While 
for most representations only one-directional translation is provided, the user 
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can change forth and back between collected form and open form, which is 
indicated in the diagram below 



. « / \ . \ . / \ . \ RemoveNCContainer 

NCContamer(UbarSpa(p), 1) * NCContamer(USpa(q), 1) — > 

UbarSpa(p) * USpa(b) SpC !^ ect Spaa(p, q) W Spa2(p, q) (23) 

SpClose 

The first in this collection of procedures is RemoveNCContainer. Its effect 
amounts to stripping off the function NCContainer and ordering the Dirac 
matrices and spinors according to the order of the spinor indices. Where nec- 
essary the flipping rules for Majorana fermions are applied (see Section |2~4"]) . 



^Procedure RemoveNCContainerQ 



Translates the indexed notation into non-commuting objects. 



At the level of non-commuting objects the dimension splitting and the 
't Hooft algebra take place, which has been implemented in the routine 
tHoof tAlgebra. 

After these steps, typically, one wants to proceed in collected form. This 
form is better suited for contractions across different spinor lines, as in 
(pi|7 /i |p 2 ](P3|7/i|P4], which can be detected very easily using argument lists 
but is very difficult to be carried out using products of non-commuting ob- 
jects. 

Traces of the form trL ■ • • trR can be converted to collected form using 
the procedure SpTrace4, all remaining products of non-commuting objects 
can be transformed by the procedure SpCollect. 



^Procedure SpTrace4(/c l5 



Operates on non- commutative objects in products representing 
a trace (trL * . . .* trR). 

Parameters 

ki , k 2 , . . . : an optional list of vectors. If 
specified, these vectors are assumed 
lightlike and traces are opened at those 
positions. If no such positions are found 
the trace is evaluated using Form's Trace4 
command. 
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The following example demonstrates the difference in treatment between 
light-like vectors and those which have not been specified as light-like. 

Local testl = trL* Sm4 ( kl ) * Sm4 ( k2 ) * Sm4 ( k3 ) * Sm4 ( k4 ) * trR ; 
Local test2 = trL*Sm4 ( kl ) * Sm4 ( k2 ) * Sm4 ( q) * Sm4 ( k4 ) * trR ; 
#call SpTrace4(q) 
Print ; 
. end 

* testl = 4*kl . k2*k3 . k4~4*kl . k3*k2 . k4+4*kl . k4*k2 . k3 ; 

* test2 = Spab (q ,k4 ,kl ,k2 , q) + Spba(q ,k4 ,kl ,k2 , q) ; 

Very often specifying the list of light-like vectors leads to more compact 
expressions, especially if the number of terms generated by taking the trace 
is very big. 

^Procedure SpCollectQ 

Translates non-commuting objects into collected form. 

This procedure requires that the n-dimensional algebra has 
already been carried out by a call to tHoof tAlgebra and that 
all massive spinors have been replaced using the routine 
LightConeDecomposition. 

As a result the all spinorial objects are expressed in terms of 
the function Spaa, Spab, Spba and Spab. 

Before one can go from collected form to open form one should eliminate 
all Lorentz contractions by the use of SpContract and SpContractLeviCivita 
The latter one is only needed if the Levi-Civita symbol e ltvpa has been intro- 
duced by taking a trace involving 75. 

The conversion from collected to open form is performed by the procedure 
SpOpen. The reverse operation, converting from open form back to collected 
form is done by the procedure SpClose. 
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^Procedure SpOpen(/ci, k 2 , • • • ) 

Translates an expression in collected form into open form. 

In open form, only the functions Spa2 and Spb2 are used to 
express spinor products. 

This procedure requires that all Lorentz indices inside spinor 
lines have been removed using the routine SpContract. 

Parameters 

ki, hi, ■ • ■ : an optional list of vectors. If 

specified, the spinor lines are only opened 
at positions indicated by the given vectors. 
In this complete translation into 

Spa2 and Spb2 functions might not 
The parameter list of SpClose allows to close spinor strings at the posi- 
tions of certain vectors, indicated by the arguments. In both cases, SpOpen 
and SpClose, the parameter lists facilitate the use of this library for unitarity 
based methods, as will be shown in the examples of Section 14.31 

^Procedure SpClose(fci, k 2 , ...) 

The inverse operation of SpOpen. Replaces the functions Spa2 
and Spb2 into the functions Spaa, Spab etc. 

Parameters 

k\ , ki, . . . : an optional list of vectors. If 

specified, only those positions are closed 
which are indicated by the given vectors. 

In some cases this operation requires that a spinor string is broken up at 
another position. Consider the following example: 

Local expression = Spab(q, kl , k2 , kl , q); 

#call SpOpenQ 
#call SpClose(q) 

Print ; 

. end 

* test =- Spab(kl , q , kl )* Spa2 (kl ,k2)*Spb2(kl ,k2); 
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However, it is not guaranteed that the operation always succeeds: if we 
had written T^call SpClose(q,kl) instead of #call SpClose(q) the following 
result would have come out, where the vector q still appears in the spinors. 

* test = Spab (k2 ,kl ,k2)* Spab (q ,kl , q) ; 



3.4-4- Contraction of Lorentz Indices 

This Section describes three routines which reduce the number of explicit 
appearances of contracted Lorentz indices. The first type of contractions, 
those involving only spinor lines, as in (pi | • • • 7 M • • • |p 2 ] (Ps\ • • • 7m ' ' ' IP4] or 
in (pi \ ■ ■ ■ 7^ • • • 7^ • • • \p 2 ], are treated by the procedure SpContract. 

^Procedure SpContract () 

Applies the equations of Section 12.31 in order to eliminate all 
Lorentz indices inside spinor lines. This procedure works on 
expressions in collected form. 

Contractions involving metric tensors are simplified by the procedure 
SpContractMetrics. This procedure is implemented such that it can deal 
with expressions in collected form but also with non-commuting objects. 

^Procedure SpContractMetricsQ 

Removes spurious appearances of the metric tensors g^, g^ u 
and g^ v . Whereas M — 4 is substituted immediately, all other 
instances of the dimensions gfi and g£ are not replaced. 

The last one in this category of procedures is SpContractLeviCivita. 
This routine uses the fact that the Levi-Civita tensor e^ upa can be written as 
a trace: 

e» vp ° = -\tr{l,Yl v l P Y} ■ (24) 

The procedure considers two cases. If the Levi-Civita tensor is contracted 
with a Dirac matrix inside a spinor string it uses the Chisholm identity to 
simplify 

= -- (n + - n_) [fT7 P - fW] (25) 

The second case is the contraction of the Levi-Civita tensor with a lightlike 
momentum. Here one can rewrite the e-tensor as 

= -~ (\p\rrr\p) - (Pirrrbi) (26) 
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Products of multiple e-tensors are reduced by applying the determinant 
relation 

^wwiw^w^ = det (gWi)*. =1 (27) 



^Procedure SpContractLeviCivita(/ci , k 2 



Eliminates Levi-Civita tensors as far as possible applying 
Equations (125)1 . (EE) and pT)). 

Parameters 

k\ , k 2> . . . : an optional list of vectors. If 
specified, these vectors are assumed 
light like. 



3.4-5. Schouten Identity 

In analogy to the package S@M [4l[ we have implemented three versions 
of the Schouten identity 

(ij)(kl) = (il)(kj) + (ik)(jl) (28) 
and its conjugated version. 



^Procedure Schouten(j9i, p 2 , qi, 52) 



Substitutes according to the Schouten identity 

{PlP2){qiq2) ->■ (Pl92)(?lP2) + (Pl?l)(P2g2)- 



^Procedure Schouten(pi, p 2 , q) 



Substitutes according to the Schouten identity 

W ■ {PiP2){qq') -> M){qT>2) + (Pi^OW)- 



^Procedure Schouten(g) 



Substitutes according to the Schouten identity 

w (qpi) _^ (p&i) (PsPx) 

Vpi,p 2 ,P3 : 1 — n — r 



(qP2)(qP3) (qP2)(P2P3) (qP3)(P2P3) 
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3.4-6. Miscellaneous Routines 

In this section we describe some routines which are of less common use. 
The procedure SpClear nullifies all expressions which should vanish due to 
spinor lines of the wrong length. The procedure SpClear is called inside 
SpContract because Equations (IT51) hold only if the spinor chains have cor- 
rect length. 

^Procedure SpClear () 

Eliminates all terms which contain a spinor line or 

where T is a product of an odd number of Dirac 
matrices, or a spinor line or [p|r|g), where T is a 

product of an even number of Dirac matrices. 

This procedure acts on expressions in collected form. 

The routine SpCheck acts in a similar way as SpClear does. Instead of 
removing the terms which should vanish trivially it marks them with the 
function SpERRORTOKEN. 

^Procedure SpCheck() 

Substitutes all terms which contain a spinor line or 

where T is a product of an odd number of Dirac 
matrices, or a spinor line (p|r|g] or [p|r|g), where T is a 
product of an even number of Dirac matrices inside the 
arguments of the function SpERRORTOKEN. 

This procedure acts on expressions in collected form. 
The last procedure of this section, SpOrder, shuffles Dirac matrices into 
a given order using the 't Hooft-Veltman algebra. The order is specified by 
the argument list of SpOrder. This routine is used in the test programs to 
bring the results into a canonical form which is necessary in order to check 
that all tests hold. 
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^Procedure Sp0rder(fci,A;2,. . . 



Sorts all products of Dirac matrices according to the ordering 
specified by the argument list. This routine acts on 
non-commuting objects. 

Parameters 

ki , k 2 , . . . : a list of vectors. If ki is to the left of 
kj in the argument list then the procedure 
shuffles fa to the left of ftj in the 
expression. 



3.5. Working with Majorana Spinors 

3.5.1. On the Relative Sign of Feynman Graphs 

The flipping rules have been tested in Golem-2.0 where the diagrams are 



generated by Qgraf [47]]. The relative sign calculated by Qgraf is incorrect 
when dealing with Majorana fermions. Here, we present a method to calcu- 
late it. Firstly we calculate (— l) p using the following code: 

Function NCOrder ; 

Id fDUMMYl?{UbarSpa , UbarSpb } (vDUMMYl?) = 

NCOrder (vDUMMYl) *fDUMMYl (vDUMMYl) ; 
Id fDUMMYl?{USpa,USpb} (vDUMMYl?) = 

fDUMMYl (vDUMMYl) * NCOrder (vDUMMYl ) ; 
#call tHooftAlgebra 
#call SpCollect 
Chainln NCOrder ; 
AntiSymmetrize NCOrder ; 
Id NCOrder (? all) = 1; 



We multiply our diagram by a non-commuting function of the external 
momenta which encodes the order of the spinors in the diagram. The argu- 
ments of this function are then brought into Form's natural ordering. The 
exchange of any two arguments results in a minus sign. 

Secondly we must determine (— 1) L with L being the number of closed 
fermion loops. This is easily done by counting the number of appearances of 
trL, as long as one ensures that the only source of spinor traces in the first 
place are closed fermion loops; the easiest way of multiplying the amplitude 
by the correct sign could be implemented as Id trL = -trL. 
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3.5.2. Fixing Fermion Chain Order 

The method we have described relies on being able to fix the fermion 
chain order. This is achieved in the procedure NCContainer. In our code an 
incoming (outgoing) Majorana fermion will initially be treated as in incoming 
(outgoing) Dirac fermion (as opposed to an anti- fermion). When we join an 
incoming Majorana spinor with an incoming Dirac fermion or an outgoing 
Dirac anti-fermion, one spinor will need to be flipped. The same applies 
when an outgoing Majorana fermion is joined to an incoming anti-fermion 
or an outgoing fermion. As an example we have an expression: 



NCContainer ( UbarSpa ( kl ) *Sm( i 1 ) * Proj Plus * UbarSpa 

>- 


(k2)) 




which is transformed using RemoveNCContainer to 


NCContainer ( UbarSpa ( kl ) , Sm(il), ProjPlus , 

SpFlip (UbarSpa( 


k2))) 




and then using (12T!) we have the fermion chain: 


UbarSpa (kl)*Sm(il )* ProjPlus *USpb(k2). 







This is the default behaviour of spinney; one can prevent spinney from 
applying the flipping rules by defining the preprocessor variable NOSPFLIP 
before calling RemoveNCContainer. 

4. Examples and Applications 

4-1. Feynman Diagram Based Reduction of One-Loop Amplitudes 

In this section we show how the Form library spinney can be used for the 
reduction of one-loop diagrams, both with conventional reduction of tensor 
integrals and with a reduction at the integrand level as described in Q. In 



the first case we use the conventions of 481. However, the method is not 



restricted to to this particular tensor decomposition. 

In our setup, we use the diagram generator QGraf 47]] . For the diagram 
in Figure [T] we obtain an output similar to the one given below 



Local diagraml = 

inp( — 1, ivlrlLl , kl , me) * 

For the sake of simplicity we skip the discussion of the color algebra. 
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Figure 1: Three-point diagram (e + e — s> tt) discussed in the example. 

inp( + l, ivlr2Ll , k2 , me) * 
out( + l, iv2rlLl , k3 , mT) * 
out(-l, iv3r2Ll , k4 , mT) * 



vertex ( 


field 


. ep ] , kl , ivlrlLl , 




field 


.em] , k2 , ivlr2Ll , 




field 


.A] , -kl-k2 , ivlr3L2) * 


vertex ( 


field 


.Ubar] , -k3 , iv2rlLl , 




field 


.U] , Q+k3 , iv2r2Ll , 




field 


.g] , -Q, iv2r3L2) * 


vertex ( 


field 


.Ubar] , -Q+k4 , iv3rlLl , 




field 


.U] , -k4 , iv3r2Ll , 




field 


.g] , Q, iv3r3L2) * 


vertex ( 


field 


.Ubar] , -Q-k3 , iv4rlLl , 




field 


.U] , Q-k4 , iv4r2Ll , 




field 


.A] , kl+k2 , iv4r3L2) * 


prop ( + 2 


-kl- 


k2, 0, iv4r3L2 , ivlr3L2) * 


prop ( + 2 


-Q, 


0, iv3r3L2 , iv2r3L2) * 


prop ( + 1 


Q+k3 


, mT, iv4rlLl , iv2r2Ll ) * 


prop ( + 1 


Q-k4, mT, iv3rlLl , iv4r2Ll); 



The first step is the substitution of the Feynman rules. In the replacement 
of the wave functions the approximation m e = is applied and helicities are 
assigned to the particles. 

Id inp ( — 1 , ivl ? ,kl ,me) = NCContainer (Ubar Spa (kl ) , ivl ) ; 
Id inp( + l,ivl ? ,k2 , me) = NCContainer ( USpb ( k2 ) ,ivl ); 
Id out ( + l,ivl ? ,k3 ,mT) = NCContainer (UbarSpb (k3 , + 1) , ivl ) ; 
Id out (-1, ivl ? ,k4 ,mT) = NCContainer (USpa(k4, -1) , ivl ) ; 
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As discussed in Section 13.2. II the letters 'a' and 'b' distinguish the helicities, 
whereas the parameter ±1 at the right-hand side differentiates between u- 
and v-spinors. 

In a similar manner the propagators and vertices are substituted. Again, 
focussing on the Lorentz algebra we omit color structure, coupling constants 
and denominators here. 

Id prop(l, kl?, mT?, ivl?, iv2?) = 

i_ * (NCContainer (Sm(kl ) , iv2 , ivl) 
+ mT * NCContainer (1 , iv2 , ivl)); 
Id prop(2, kl?, 0, ivl?, iv2?) =- i_ * d(ivl, iv2 ) ; 
Id vertex ([field .ep]? , kl?, ivlLl?, 
[field. em]?, k2?, iv2Ll?, 
[ field .A]? , k3? , iv3L2?) = 
NCContainer (Sm(iv3L2 ) , ivlLl , iv2Ll): 



Now, a call to RemoveNCContainer brings the non-commuting object in the 
right order and simplifies the representation of the expression. We replace 
the massive spinors by projecting onto massless vectors. Hereby we introduce 
new vectors I3 and U as defined in Equations (0) and (jSJ). 

N 

T^call RemoveNCContainer 

T^call Light ConeDe compos it ion (k3 , 13 , kl ,mT) 
T^call LightConeDecomposition (k4 , 14 ,k2,mT) 



At this point the manipulation of the diagram branches depending on 
the output one wants to achieve. In a conventional reduction of the tensor 
integrals the algorithm contiues after replacing the integration momentum 
Q by the corresponding expressions for the tensor integrals. In the case of 
a reduction at the integrand level we like to express Q in terms of a four- 
dimensional projection Q4 and a scale /1 2 such that Q 2 = (Q4) 2 — /i 2 . 

4-1.1. Conventional Tensor Reduction 

The expression corresponding to the Feynman diagram in Figure [T] has 



integrals of rank two at most. Following the notation of |48[ the momenta in 



the loop are r\ = —k.4, r2 = k% and r% = 0. The corresponding form factor 
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representation of the occurring tensor integrals is 

JJ(5) = A 3 '°(S); (29) 
I^(S) = -k^Al'^S) + k^A 3 /(S); (30) 

i 3 l^DJ — K 4 K 4 I^DJ ^/C 3 K 4 -+- K 4 K 3 J/i 12 I^DJ 

+ fcfC4 2 (5) + ^ M1M2 5 3 ' 2 (S). (31) 

These formulae are easily implemented for the given example. It should 
be noted that the calls to the function LightConeDecomposition must be 
repeated after this step since we reintroduce the momenta k$ and k± here. 

N 

ToTensor, Functions, Q, QTens ; 

If (count (QTens,l) = = 0) Multiply A30 ; 

Id QTens(ivl?) = -k4(ivl) * A31(l) + k3(ivl) * A31(2); 
Id QTens(ivl?,iv2?) = 

+ k3(ivl)*k3(iv2) * A32 ( 1 , 1 ) 

- k3(ivl)*k4(iv2) * A32(l,2) 

- k4(ivl )*k3(iv2) * A32(l,2) 
+ k4(ivl )*k4(iv2) * A32(2,2) 
+ d(ivl , iv2) * B32 () ; 



The remaining steps of the program carry out the 't Hooft algebra and 
simplify the expression by taking out all contractions of Lorentz indices. 
The replacements of dEps are consistent with the dimension splitting in use 
and define the symbol eps such that n = 4 — 2s. In the four-dimensional 



helicity (FDH) scheme |49H53l| the symbol eps can be replaced by zero; in 
the 't Hooft- Veltman scheme the products between s (resp. s 2 ) and the form 
factors, which are formally Laurent series in s, lead to rational terms. After 
these steps the diagram is in a form suitable for numerical evaluation if one 
provides routines for the computation of spinor brackets ((■■) and [■■]) and a 



library of integral form factors such as golem95 54 



#call tHoof tAlgebra 

#call SpCollect 

T^call SpContract 

T^call SpContractMetrics 

Id dEps(Q?{kl ,k2 , 13 ,14}, ivl?) = 

Id dEps (ivl?, ivl?) = -2*eps ; 

T^call SpOpen 
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4-1-2. Reduction at the Integrand Level 

Following the strategy of reducing the Feynman diagrams at the integrand 
level, we show in this section how the numerator of the Feynman diagram can 
be constructed very easily by the use of spinney. We consider two different 
approaches corresponding to two publicly available reduction packages. The 



strategy implemented in CutTools [55] uses a numerator function N(Q 4 ) 



which only depends on the four- dimensional projection of the (complex) in- 
tegration momentum Q. This method allows for the reconstruction of the 
cut-constructible terms but only partially recovers the rational parts of an 
amplitude. This is due to the fact that terms in fi 2 , where Q 2 = Ql — fi 2 , lead 
to rational terms which are not taken into account and need to be added by 



a separate calculation [56|, |57| . An improved reduction at the integrand level 



has been implemented by the authors of Samurai [15|. This method takes 
advantage of a numerator function depending on both Q 4 and /j, 2 , which can 
be decomposed into 

A(Q 4 , /i 2 ) = N Q (Q 4 , /i 2 ) + eJVi(Q 4 , /i 2 ) + e 2 N 2 (Q 4 , /x 2 ) + 0(e 3 ). (32) 

As in the case of the conventional tensor reduction, e is set to zero (keeping 
terms in /1 2 ) in the FDH scheme, whereas in the 't-Hooft Veltman scheme 
the independent reduction of the numerators N , Ni and N 2 , multiplication 
with the apropriate terms of the Laurent series of the scalar integrals leads 
to the full result, including both both the cut-constructible and the rational 
part of the amplitude. 

Irrespective of the approach chosen the first step to be taken should be 
to carry out the 't Hooft algebra. Then the two methods differ in the way 
they deal with terms in dEps. Using CutTools these terms are neglected 
and set to zero. In the case of Samurai one substitutes the relations which 
are implied by Q 2 = Q\ — fi 2 , where g^Q 2 = Q^Q V and Q 4 = g^Qu, and 
therefore fi 2 = —g^Q^Qu- 

#call tHoof tAlgebra 
#call SpCollect 
#call SpContractMetrics 
Id dEps(Q,Q) = - mu2; 
Id dEps(Q?, ivl?) = 0; 
Id dEpsjivl?, ivl?) = -2*eps ; 
Argument Spaa , Spab , Spba , Spbb , d4 ; 
Id Q = Q4; 
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End Argument ; 

Id Q.Q = Q4.Q4 - mu2; 
Id Q = Q4; 

Id d4(kl?, ivl?) = kl(ivl ); 

Finally, one can contract the remaining Lorentz indices and bring the the 
expression into a form where spinor brackets ((••), [••] and [-| Qi\-)) are the 
only functions to be evaluated numerically. 

T^call SpContract 

#call SpOpen(kl ,k2 , 13 , 14 ) 



4-2. Working with Majorana Spinors 

4-2.1. Charge Conjugation of a Vector Current 

We illustrate our approach through the following simple example. The 
vector current, (p + \ r y' J '\q + ) , satisfies the charge conjugation relation 

(p+lfV) = (q-\r\p-), (33) 

which we can show explicitly using charge conjugation relations. This charge 
conjugation operation is equivalent to us reversing the fermion flow arrow. 
Therefore we can use Equations (I20I) to show this equivalence holds: 

(p + \r\q + ) = W\q) {-l) P (q\r'\v\ = ("1) P (-1) (q\l>} (34) 
= (q\l>] = (q-\l>-) (35) 

where we have used (|T9l) and (12 ID and we have defined our reference order as 
(p,q) giving (— l) p = —1. This result is true for both Dirac and Majorana 
fermions. 

4-2.2. Majorana Exchange 

We now consider two Dirac fermions scattering to two vector bosons via a 
t-channel Majorana fermion exchange. There are two equivalent orientations 
to choose from as shown in Figure [2j We write our amplitude as 

A=A»%(p 3 )e;(p 4 ). (36) 

Our first choice of orientation gives 

Ar=( P i\WS(p)Y\P2) (37) 
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Figure 2: The process (// — > VV) mediated by the exchange of a Majorana fermion (x) 
as discussed in the example. The left hand figure shows the original diagram; the middle 
figure shows the first choice of fermion orientation; and the figure on the right shows the 
second choice of orientation. 



and the second orientation gives 

AT = (-1)(P2\(Y)'S(-p)Y\Pi) (38) 

where we have chosen the reference order (p\,p2). Applying our flipping rules 
to Equation (1371) gives the result: 

Af = (-l) P (P2\(Y)'S(- P )r\Pi) = A?. (39) 

Therefore the amplitude is independent of the original choice of fermion 
orientation. 

In our code, part of the output is: 

NCContainer (USpa( kl ) * SpFlip (Sm(mu) * Proj Plus ) *Sm( k4 ) 

*Sm(nu)*ProjPlus*USpb(k2))*inv(es23 ) . 
> 

Upon applying the RemoveNCContainer routine we obtain the result 

N 

UbarSpb (kl )* ProjPlus*Sm(mu)*Sm(k4)*Sm(nu) 

* Proj Minus *USpa(k2 ) * inv ( es23 ) 
> 

We have picked up a minus sign from the flipping of the 7 M . What remains 
is to multiply by (— l) p as explained previously 

4-3. Coefficients of Scalar Integrals by Unitarity Based Methods 

In this example we consider Bhaba scattering in QED at the one-loop 
level. We use the well known fact that any (leg-ordered) one-loop amplitude 
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can be written in terms of a basis of scalar one-loop integrals, 



D~ + % Kn J D~D~ + . S . c --3 J 



d n l 



Di 2 Di 3 



ll H<l2 " U<«2<«3 

[ d n l 

+ 2^ d W3uJ D D D D , (40) 

il<i 2 <i3<i4, 1 2 3 4 



where A = [(/ + r^) 2 — m 2 ]. 

We consider the two diagrams in Fig. [3j The box coefficient is isolated by 
performing four cuts, which completely disconnect the loop amplitude into 
4 tree- level partial amplitudes with three external legs each. In fact, since 
there is only one vertex in QED and the identity of one of the legs is fixed (it 
is one of the original external fermions), we have two possible contributions 
in each amplitude differing by the exchange of the remaining fermion / photon 
legs of the QED vertex. Consistency, i.e. the fact that the lines connecting 
neighbouring vertices need obviously be of the same type, reduces this to two 
choices, depicted on the left-hand side of Fig. |3j The helicity conservation 
on the fermion lines and the restrictions of complex kinematics allow us to 
write down the only four possible internal helicity configurations. The four 
cuts also provide four independent constraints on the loop momentum, thus 
determining all its components (in fact there are two solutions). For details 



sec 



22|. We take all momenta to be outgoing. 
We begin the computation by preparing a small procedure that will be 
used repeatedly together with SpOpen to simplify denominator structures. It 
is an ad-hoc solution based on the knowledge of type and depth of denomi- 
nators arising due to loop momenta substitutions. 

N 

^Include— spinney, hh 
^Procedure SpSimplifyQ 
Repeat ; 

FactArg SpDenominator ; 

ChainOut SpDenominator; 

Id SpDenominator (SpDenominator (Spa2 (vDl? ,vD2?))) = 

Spa2(vDl,vD2); 
Id SpDenominator (SpDenominator (Spb2 (vDl? ,vD2?))) = 

Spb2(vDl,vD2); 
Id SpDenominator ( SpDenominator (fDl? ( ? vDl) ) ) = 

fDl(?vDl): 
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Id SpDenominator (sDUMMYl? number _ ) = 1/sDUMMYl; 
EndRepeat ; 
#EndProcedure 

> ; 

We write down the explicit formula for one of the four internal helicity con- 
figurations, split into separate numerator and denominator expressions. The 
formulae for the remaining configurations are similar. 

Global FaNUM = -4 *Spbb(ll , 10) *Spaa(10,13) *Spaa(ll,12) 
*Spbb(13,12) *Spaa(pl,R0) *Spbb(R0,p2) 
*Spaa(p3,R2) *Spbb(R2,p4) ; 

Global FaDEN = Spaa(ll,R0) *Spbb(R0,ll) *Spaa(13,R2) 
*Spbb(R2, 13 ) ; 

The first step is to substitute the solutions for the loop momentum on the 
cut, which in our case is: 

* = (41) 
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In order to do so we need to re-express all loop momenta in terms of l and 
external momenta. We use SpClose to have the spinor products in a form 
which is suitable for standard linear substitutions: (p\ . . . 1^ . . . k 2 . . . \q) , with 
p and q being external momenta and U the loop momentum. This might not 
be immediately possible, as in our example, where SpClose produces an 
output of the type (ZjJ . . . \l i2 ) instead. However, a call to Schouten solves 
the problem at the cost of having more terms in the calculation. 



#Call SpOpen 




#Call Schouten(10 , 13 ,pl ,R0) 




#Call SpClose (10 ,11 ,12 , 13)' 




. sort 




ToTensor Functions 11 , LI; 




ToTensor Functions 12 , L2 ; 




ToTensor Functions 13 , L3 ; 




Id Ll(iDl?) = 10 (iDl) - pi 


;iDl); 


Id L2(iDl?) = 10 (iDl ) - pi 


;iDl) - P 2(iDl); 


Id L3(iDl?) = 10 (iDl) + p4 


;iDl); 


#Call SpContract 




. sort 





Now we substitute the solution for /q and contract all explicit Lorentz indices. 

ToTensor Functions 10 , L0 ; 
ChainOut L0 ; 
. sort 

Id L0(iDl?) = Spab(pl ,iDl , P 4)*Spb2(p2 ,pl)*(l/2)* 

SpDenominator ( Spb2 (p2 , p4 ) ) ; 
#Call SpContractMetrics 
#Call SpContract 
#Call SpOpen 
. sort 

We combine the numerator and denominator expressions obtained before. 
Our formulae still involve the reference momenta of the photons, which is 
our gauge freedom - we can greatly simplify the computations by making an 
explicit choice at this point. 

Global Fa = FaNUM* SpDenominator (FaDEN) ; 
Argument ; 
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Id R2 = pi ; 
Argument ; 

Id R2 = pi; 
EndArgument ; 
End Argument ; 
#Call SpOpen 
#Call SpSimplify 
#Call SpOpen 
. sort 

As a finishing touch, we apply the following sequence of transformations to 
get rid of any remaining double denominator-type expressions, which appear 
upon introducing the explicit form of the Iq loop momentum. 

Argument ; 

#Call Schouten(pl,p2,p4,R0) 
EndArgument ; 
Repeat ; 

#Call SpOpen() 

Argument ; 

#Call SpOpen () 

EndArgument ; 

#Call SpSimplify 
EndRepeat ; 
. store 

Finally, we write down the expression for the complete box coefficient. 
Note the factor 1/2, which comes from averaging over two solutions to the 
loop momentum constraints. 

Local F = (1/2)* (Fa + Fb + Fc + Fd); 
Print F; 
. end 



The FORM output we obtain is: 

, = 2(p 1 p 2 ) 3 (p 3 p 4 )[p 2 p 1 ] 2 2(p l p A ) 3 (p 2 p 3 )[p 4 p 1 } 2 
(P2P4} 2 (P2P4} 2 

+ 2(pip 2 ) (piP3> \p2P1] [P4P2] - 2(pip 3 ) (pm) \p4p1} [pm]- (42) 
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5. Conclusion 



In this article we have presented a new Form library for processing expres- 
sions containing helicity spinors in four and n = (4 — 2e) dimensions. For 
the n-dimensional algebra we have implemented the 't-Hooft algebra with 
dimension splitting. This gives full flexibility to the user about the choice of 
the regularisation scheme, as in many schemes different from the 't Hooft- 
Veltman scheme it is sufficient to neglect terms in e which stem from the 
numerator algebra. 

In various examples we have shown that the new package is applicable 
to calculations both using traditional and modern, unitarity based methods. 
The provided routines implement the typical steps which are necessary for 
an algebraic simplification of helicity amplitudes and therefore simplify the 
task of implementing such calculations. The implementation of flipping rules 
for Majorana spinors allows to extend the domain of applicability to theories 
beyond the Standard Model. 
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